Coexistance of giant tunneling electroresistance and magnetoresistance in an all-oxide 

magnetic tunnel junction 
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We demonstrate with first-principles electron transport calculations that large tunneling magne- 
toresistance (TMR) and tunneling electroresistance (TER) effects can coexist in an aU-oxide device. 
The TMR originates from the symmetry-driven spin fihering provided by the insulating BaTiOa bar- 
rier to the electrons injected from SrRuOa. In contrast the TER is possible only when a thin SrTiOa 
layer is intercalated at one of the SrRuOa/BaTiOa interfaces. As the complex band-structure of 
SrTiOa has the same symmetry than that of BaTiOa, the inclusion of such an intercalated layer 
does not negatively alter the TMR and in fact increases it. Crucially, the magnitude of the TER 
also scales with the thickness of the SrTiOa layer. The SrTiOa thickness becomes then a single 
control parameter for both the TMR and the TER effect. This protocol offers a practical way to 
the fabrication of four-state memory cells. 



Epitaxial magnetic tunnel junctions (MTJs), display- 
ing giant tunneling magnetoresistance (TMR) at room 
temperature [H Q , represent the enabling technology for 
ultra-high density magnetic data storage. In MTJs the 
insulating barrier plays a dual role; it magnetically de- 
couples the two ferromagnetic (FM) electrodes so that 
their magnetizations, M, can be arranged either paral- 
lel or antiparallel to each other, but it can also act as a 
spin filter, if epitaxially grown. This is due to the wave- 
function symmetry selective decay of tunneling electrons 
across a crystalline insulator. As the two spin-manifolds 
of the Fermi surface of a FM metal present different sym- 
metries, such symmetry selectivity translates into spin 
selectivity. This is the case for Fe/MgO(001), where the 
barrier is more transparent to the tunneling of electrons 
with Ai symmetry. These are present in Fe only for the 
majority spin 0,131, so that the Fe/MgO(001) stack ef- 
fectively behaves as a half-metal. 

In conventional MTJs, however, the insulating barrier 
is a passive element, i.e., its electronic structure can- 
not be changed by external stimuli. A different situa- 
tion is encountered when using a ferroelectric (FE) ma- 
terial. A FE is intrinsically insulating and at the same 
time possesses a macroscopic order, the electrical polar- 
ization, P. When FE materials are incorporated into a 
tunnel junction, one expects the junction resistance to 
become dependent on the direction of P with respect to 
the layer stack, an effect known as tunneling electrore- 
sistance (TER) 0. It then becomes natural to think 
about devices combining materials with both FM and 
FE orderings Q. Here, one can exploit the possibility 
of manipulating the two independent order parameters, 
P and M, by means of their conjugate fields, namely 
the electric and magnetic fields. The fabrication of FE 
random access memories with non-destructive reading 
demonstrates the potential of such an approach. 

Although it is possible, at least in principle, to ob- 
tain a large TMR in MTJs with FE barriers [1, Q, it is 
sensibly more complicated to obtain a large TER. The 



key ingredient for a MTJ to show TER is that it should 
exhibit inversion symmetry breaking. This is almost al- 
ways the case in real devices as unintentional disorder 
breaks the symmetry. However, disorder is hardly con- 
trollable. Furthermore, even when the entire junction 
termination is different at either side of the insulating 
barrier the TER appears relatively modest Q. A sec- 
ond strategy uses two magnetic electrodes made of dif- 
ferent metals and thus different abilities to screen surface 
charges Q. In typical metals with high carrier mobility, 
however, the screening length is short and the surface 
charge is strongly localized at the interface. The result- 
ing potential profile thus remains approximately mirror- 
symmetric upon polarization reversal and the expected 
TER is small. Even for Fe/BaTiOa/Lao.eTSro.gsMnOa 
junctions, where Lao.gTSro.asMnOa is quite a poor metal, 
a TER of only 37% has been observed 0. 

One can bring the concept of having different screening 
lengths at the two sides of the FE layer to the extreme 
by including a second insulator (INS) in the stack, i.e., 
by considering a FM/INS/FE/FM junction. An interest- 
ing example of such structure is when the INS is actually 
vacuum, as in a scanning tunnel microscopy experiment. 
Here, an extremely large TER has already been measured 



ll|,|l2|. In this letter we demonstrate, by first principles 



electronic transport calculations, that the same effect can 
be achieved in an all-oxide solid state device. In par- 
ticular, we show that intercalating a few monolayers of 
SrTiOs into a SrRuOs/BaTiOs/SrRuOs junction creates 
an additional potential barrier that is switchable with 
the FE polarization of BaTiOa. This not only amplifies 
the TER, but also makes it exponentially dependent on 
the SrTiOa thickness. Furthermore, as SrTiOa is elec- 
tronically very similar to BaTiOa and thereby provides 
comparable spin-filtering for SrRuOa [1,0], the junction 
also displays a remarkably large TMR. 

The electronic structure of the junction is calculated 
by using density functional theory (DFT) as numerically 
implemented in the siesta code [l3[. Structural relax- 



ation is performed with the Perdew-Burke-Ernzerhof gen- 
erahzed gradient approximation (GGA) to the ex- 
change and correlation fmictional, with a 6x6x1 fc-point 
Monkhorst-Pack mesh and a grid spacing equivalent to 
a plane- wave cutoff of 800 eV. In contrast, for the elec- 
tronic properties and the transport we use the atomic 
self-interaction correction (ASIC) scheme built over 
the local spin density approximation. ASIC has been pre- 
viously found to improve the electronic properties of bulk 
BaTiOa [iF 



and SrRuOa 



17| and it is vital in transport 
calculations where one has to ensure a good band align- 
ment between dissimilar materials [18|. Unfortunately, 



even in its variational form [19| the ASIC scheme does not 
describe the FE phase accurately enough, so that a com- 
promise is required; the GGA is used for the relaxation 
and the ASIC for the transport calculations. Electron 
transport is computed with the ab initio code SMEAGOL 



-22|, which combines DFT with the non-equilibrium 
Green's functions scheme. SMEAGOL uses siesta as its 
DFT platform so that the same convergence parameters 
are employed for the transport, except for the fc-point 
sampling where we consider a much larger 100x100x1 
mesh. 

The supercell considered here comprises 6 BaTiOa 
unit cells 2.5 nm) and 3 unit cells of SrRuOa at- 
tached at each side to function as electrodes. Fur- 
thermore, we intercalate a thin SrTiOs layer between 
BaTiOa and SrRuOa at one side of the junction so that 
that the final stack is (SrO-Ru02)3/(SrO-Ti02)™/(BaO- 
Ti02)6/(SrO-Ru02)3, where m = 0,1,2. The in-plane 
lattice parameter is set to that of bulk SrTiOa (3.95A) to 
mimic the effect of a SrTiOs substrate. This applies com- 
pressive strain to both SrRuOa and BaTiOa and in doing 
so increases the polarization of BaTiOa. The Berry phase 
method gives a GGA polarization of 43.8//C/cm^ for the 
bulk (c/a = 1.05) and 48.1yU C/cm^ for the strained struc- 
ture (c/a — 1.08). Note that the GGA systematically 
overestimates the polarization of FE oxides, but such a 
detail docs not affect our conclusions. 

We consider two alternative geometries for the junc- 
tion, characterized by the BaTiOa polarization pointing 
in opposite directions. Both geometries are relaxed to 
a tolerance of 40 meV/A (less than 4 meV/A for the 
m = case). When BaTiOs is included in the capacitor 
structure the displacements at the center of the supercell 
correspond to a polarization of 35.5/iC/cm^, i.e., sensibly 
reduced from its bulk value. 

An indication of the polarization structure is obtained 
from Fig. [1] where we show the atomic displacements, 
S, along the MTJ stack. Here S = (zcation — -^o), where 
2^cation and zq denote respectively the cation and O po- 
sition in a particular plane. As such, S > defines a 
structure with the polarization pointing parallel to the 
substrate normal and away from the intercalated SrTiOs 
layer (P-j.), while 5 < define a structure with the polar- 
ization pointing in the opposite direction (/V). Clearly, 
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FIG. 1. Relative cation-oxygen displacements along the z-axis 
(the junction stack direction) for the fully relaxed structure: 
(a) m = 0, (b) m = 1 and (c) m = 2. The black (red) 
line corresponds to displacements in the BO2 (AO) planes, of 
ABO3. The sohd (dashed) lines indicate P-> (P<-). 



as far as the displacement is concerned there are no 
significant differences between the intercalated SrTiOs 
and BaTiOs, which means that SrTiOa, an incipient FE, 
takes on the FE distortion of BaTiOa. This is valid only 
for thin SrTiOa films (m = 1,2) while we expect that 
thicker layers will lose the FE state. 




FIG. 2. Complex and real band structure of bulk SrTiOs 
(left panel) and bulk BaTiOa (right panel) , calculated for the 
FE structure constrained to the in-plane lattice parameter of 
SrTiOs. 

The electronic structures of SrTiOa and BaTiOa are 
also rather similar to each other as shown in Fig. [21 where 
one can observe that the real band structures of the two 
materials almost coincide. Furthermore, and more im- 
portantly here, the symmetry of the complex part of 
the band structure is identical in both materials, with 
a Ai symmetry band dominating the lower energy part 
of the band gap and a A5 one defining the region near 
to the conduction band. We then expect that interca- 
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lating a SrTiOa layer will give a MTJ with the same 
spin- filtering properties of the SrRuOs/BaTiOa/SrRuOs 
stack Q. Here, we use the "optimistic" TMR ratio, 
^TMR — '^tT^-G'ti ^ where G-^^ {G^_i) is the total con- 
ductance for the parallel (antiparallcl) orientation of the 
magnetization. The TMR is found to increase with 
SrTiOs thickness due to the increasing length of the spin- 
filtering barrier. In particular, for m = 2 at zero-bias 
-Rtmr exceeds 10^% (the actual value depending on the 
polarization direction) , meaning that at these thicknesses 
the barrier acts as an almost perfect spin-filter. 

We now discuss the TER effect in the junction by first 
looking at the electrostatic potential profile. In order to 
sustain the internal electric field associated with a FE 
material, the electrostatic potential profile must display 
a finite slope. Concurrently, assuming that no bias is 
applied across the junction and that the two electrodes 
are made of identical materials, i.e., they have the same 
Fermi energy, the average potential in the electrodes 
should be identical. As a consequence, it is necessary 
that surface charge forms at the interface between the 
FE layer and the metallic electrodes. This creates a de- 
polarizing field so that the potential across the interface 
can be matched and also sets the critical thickness for 
the onset of the FE state in a thin film (isj . 
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FIG. 3. Difference in charge distribution and electrostatic 
potential profile through the m = 2 junction. These are 
the differences between the relative quantities (charge den- 
sity and electrostatic potential) as calculated for the FE and 
centrosymmetric structures and averaged over the plane per- 
pendicular to the junction stack. The left (right) panel is for 

In Fig. [3] we present both the charge density and the 
electrostatic potential profile across the m = 2 junction. 
These are obtained as the planar average of the differ- 
ence between the relevant quantity calculated for the cen- 
trosymmetric and FE configurations. The atomic oscilla- 
tions thereby cancel and one is left with the modifications 
of the charge density and the potential due to the onset of 
the FE phase. In general we observe that charge density 
of opposite sign forms at either side of the FE layer result- 
ing in the expected potential difference. As we move into 
the metallic layers at the BaTiOa/SrRuOs interface an 
additional peak in the charge density can be seen, which 



acts as a depolarization charge and brings the potential 
back to zero. In contrast, at the SrTiOs/BaTiOa in- 
terface there are not sufficient screening charges so that 
the depolarization charge forms instead at the metallic 
SrRuOs electrode. This leaves the potential in SrTiOa 
pinned to that at the interface with BaTiOs. Thus, when 
one reverses P the potential in the SrTiOa layer is rigidly 
shifted. Note that, as a consequence of such charge dis- 
tribution, the average electrostatic potential in SrTiOs 
remains flat despite the ionic displacement. 
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FIG. 4. Total density of states (DOS) for the m = 2 junction 
(green shaded area) and DOS projected over the SrTiOs layer 
(solid black area) for the two different polarization directions, 
namely (a) P-^ and (b) P<_ . Note the rigid shift in the SrTiOs 
DOS as P is reversed. In the insets we show a schematic of 
the barrier profile for both polarization directions. 

This rigid shift in the potential can be appreciated by 
looking at Fig.|4j where we show that the density of states 
(DOS) projected onto the SrTiOa layer is rigidly dis- 
placed by the reversal of the BaTiOa polarization direc- 
tion. In particular, for the P<_ configuration the SrTiOs 
conduction band edge is considerably closer to the junc- 
tion Fermi level, Ef, than for the P-^ case. This means 
that the height of the SrTiOa potential barrier presented 
to the tunneling electrons changes according to the di- 
rection of P. In summary, the overall scattering poten- 
tial appears as follows: for there is a high barrier 
in SrTiOs followed by a triangular barrier in BaTiOs, 
which decreases as one moves away from the Ti02/BaO 
interface [see insets of Fig. IH^a)]. In contrast, for P<_ 
the SrTiOs barrier is small while the triangular BaTiOs 
barrier increases away from the interface. As a conse- 
quence, tunneling through BaTiOa is essentially insensi- 
tive to the polarization direction, as the two triangular 
barriers are identical. However, the barrier height across 
SrTiOa changes significantly with the P direction. Such 
a polarization-dependent change in the SrTiOa barrier is 
the cause of the TER effect in this junction. 

In Tabic [J we summarize our transport results. In par- 
ticular, we present the junction conductance at zero bias 
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TABLE I. Layer conductance (in units of f2 ^cm ^) and both 
TMR and TER ratios (in %) for different m. Here G^"' is tlie 
layer conductance for the magnetic configuration ocr' and the 
electrical polarization pointing in the oi direction. Note that 
the TER depends on the magnetic configuration of the junc- 
tion and the TMR depends on the electrical configuration. 



for the two different P directions (cither ■(— or — )■) and the 
two different magnetic arrangements of the electrodes (ei- 
ther parallel tt or antiparallel t^|. orientation) for m = 0, 1 
and 2. In the table, in addition to iJxMR, we also report 
the figure of merit for the TER effect, namely the TER 



ratio i?TER 



Note that the TMR is now 



dependent on the polarization direction so the ratio is 
defined as i?TMR = °^t^ ° ; being C^"^ the junction 
conductance for the aa' magnetic configuration and P 
pointing in the a direction. 

Firstly, it can be observed that our MTJ sustains a 
very robust TMR regardless of the direction of the po- 
larization vector. This is simply a consequence of the 
spin-filtering effect and of the fact that the electronic 
structure of SrTiOa and BaTiOa are rather similar. More 
interesting is the dependence of the TER on the SrTiOa 
barrier. Since in our junction the TER originates from 
a change in the SrTiOa barrier height, the effect is ex- 
pected to be magnified by increasing the barrier width, 
i.e., the SrTiOa layer thickness. This is indeed the case, 
as demonstrated by the dependence of i?TER on m re- 
ported in Table m In particular, we find that i?TER ^ 
for m = 0, i.e., when there is no intercalated SrTiOa. 
It then increases drastically for m — 1 and m = 2. 
This increase is, in fact, exponential in m, and it goes 
as e-('^"-'^")™, where A" is the SrTiOa barrier height 
in the Pa configuration. This is an important result, as 
it demonstrates that the TER can be tuned to a great 
degree by simply controlling the SrTiOs layer thickness. 

It is also important to note that for a given (m ^ 0) 
junction there are four very distinct conductance states 
depending on both the magnetization direction of the 
electrodes and the polarization direction of the ferroelec- 
tric layer. This means that our proposed device can op- 
erate as a four state memory cell with four well-separated 
conductive states. Finally, one can quantify the depen- 
dence of the TMR on the P direction by calculating the 
tunneling electro-magneto resistance (TEMR) ratio, de- 
fined as i?TEMR = ^™r~J^tmh . For m = 2 we find 

-"'TMR 

-Rtemr = 460%, a value which is comparable to those re- 



ported experimentally for Fe/BaTiOs/Lao.erSro.saMnOa 
thin film structures (ranging between 140% and 450%) 

In conclusion, we have discussed the effects of includ- 
ing a wide-band gap insulator into a MTJ based on a FE 
barrier. We have demonstrated that in such a junction 
the tunneling barrier profile can be tuned by reversing 
the direction of the macroscopic electrical polarization. 
This results in a tunable TER effect which may coexist 
with a TMR effect. In particular, the choice of SrTiOs 
and BaTiOs, which both offer excellent spin- filtering to 
spins injected from SrRuOa, results in a device which 
also displays remarkably large TMR ratios. Importantly, 
both the TMR and the TER are tunable and increase 
with the SrTiOa layer thickness. As such our proposed 
stack offers a robust protocol for constructing devices dis- 
playing simultaneous TER and TMR effects. These can 
operate as a four state memory element for data storage 
applications. 
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